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ABSTRACT 

We study the orientation and density profiles of the cosmological voids with SDSSIO 
data. Using voids to test Alcock-Paczynski effect has been proposed and tested in 
both simulations and actual SDSS data. Previous observations imply that there exist 
an empirical stretching factor which plays an important role in the voids’ orientation. 
Simulations indicate that this empirical stretching factor is caused by the void galax¬ 
ies’ peculiar velocities. Recently Hamaus et al. found that voids’ density profiles are 
universal and their average velocities satisfy linear theory very well. In this article 
we first confirm that the stretching effect exists using independent analysis. We then 
apply the universal density profile to measure the cosmological parameters. We find 
that the void density profile can be a tool to measure the cosmological parameters. 
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1 INTRODUCTION 

Observations have shown that the universe is isotropic and 
homogeneous at large scales. At the same time, the uni¬ 
verse is highly anisotropic at smaller scales and builds up 
a hierarchical structure of matter from galaxies to super 
clusters. While these dense structures are forming, it is 
unavoidable to create underdense regions which are called 
voids. Like the overdense regions, the voids are also highly 
affected by the evolution of the universe’s energy density 
compone nts and therefore could be a powerful cosmologi¬ 
cal probe dThompson fc Gregory! I 2 OI 1 I I. For example, their 
shapes and sizes are sens i tive to the nature of the dark 
energy (|Biswa ^_eL^ I l2010l : IBos et all I 2 OI 2 I : Ijennings et al.l 
12013: iLi et al.ll2012r) and their internal dynamics may repre¬ 
sent t he fifth force in the modified gravity t heory (ICai et al.l 
l2015l : [ciampitt et al.ll20l3 : lLi &: Zhadl200^ 1. 

Unlike luminous oyerdense regions, the voids can not 
be found directly, partly because they are generally faint 
(because of the lack of galaxies) or even completely dark. 
Thus, they are identified through mappi ng of the overdense 
areas first and then by using void finders dSutter et al.ir2015l l 
to pick up the voids. At large scales, the voids’ shapes are 
statically isotropic and therefore stacking could make the 
voids look like standard spheres. This is the basic require¬ 
ment to precede the Alcock-Paczynski test 
I 1979 I I. and several such tests have been done 


(lAlcoc^e£_^ 

(Sutter et al. 
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l2012bl . 2014^Q}_ei^er through the Millenium Run simula¬ 
tion i Springel et al.lBoOSl l or actual observations. 

Since voids’ shapes have been affected by the evolu¬ 
tion of the cosmological components, their orientation will 
also be affected by the evolution. Thus, we could also use 
orientation of the voids to test a gravitational model in 
question. Simulations have shown that the orientations of 
voids ar e highly correlated with their close ne i ghbors (within 
30Mp c)( Platen_et^ 20081: Lee^^ra 2000l : IPorciani et al.l 
l2002f : iLee fc Parkll2006l : IPark fc Leell2007l L This, however, 
has not been confirmed in any actual observation. In this 
article we are not trying to test this short distance align¬ 
ment, instead we are focusing on the large scale orientation. 
At such scales Alcock-Paczynski effect and galaxies peculiar 
velocities are the main source of distortion. Our study con¬ 
firms that empirical stretching factor (which is caused by 
peculiar velocities) also plays an important role in the ori¬ 
entation of the voids and their density profiles. As a conse¬ 
quence, by studying the void density profile, we could study 
the cosmological energy density components. In the follow¬ 
ing section we first desc ribe our dataset fr o m the public 
Cosmic Void Catalogue llSutter et al.l l2012al . l2014al lai and 
then show the distribution of the orientation of the void. 
The Alcock-Paczynski effect strongly affects the statistics of 
the long axis and could be a good probe to test the cosmi- 
cal paramete rs. The orientation test is similar to the one by 
Sutter et al. dSutter et al.ll2014d i. though the method is dif¬ 
ferent. However, this is not enough to quantify the peculiar 
velocity effect. Therefore, we fit the void density with an 
empirical universal density. This in turn could quantify the 
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Table 1. The void number in different section 


Sample Name 

^min 

^max 

redshift 

-^void 

comoving 

-^void 

dr72diml 

0.0 

0.05 

110 

102 

dr72dim2 

0.05 

0.1 

186 

184 

dr72brightl 

0.1 

0.15 

186 

188 

dr72bright2 

0.15 

0.2 

112 

96 

drl01owz2 

0.2 

0.3 

138 

136 

drlOlowzS 

0.3 

0.4 

215 

198 

drl01owz4 

0.4 

0.45 

71 

90 

drlOcmassl 

0.45 

0.5 

227 

229 

drl0cmass2 

0.5 

0.6 

694 

696 

drlOcmassS 

0.6 

0.7 

386 

X 


peculiar velocity effect. From the likelihood analysis 

change with SIa, we find that the density profile can serve 
as the cosmological probe, though the current void number 
might not be fully sufficient yet for very precise results and a 
prior function must be known to study the confidence level. 


2 DATASET 

Recentl y, several void catalogues were completed an d 
released l|Nadathur fc Hotchkiss! 20141 : ISutter et al.l [ 2 (fl 2 ah . 
In this st udy, we use the voi d s in th e Public Cosmic Void 
catalogue llSutter et al.l[2012al . I^Oldal Hl. because this cata¬ 
logue includes the most recently results. This Catalogue puts 
the voids and void galaxies in the following coordinates 


Xl 

= Dc( 2 ) cos(Dec) cos(Ra) 


X2 

= Dc( 2 ) cos(Dec) sin(Ra) 

( 1 ) 

X3 

= Dc( 2 ) sin(Dec) 



Here 2 is redshift. Dc{z) is the comoving distance from 
the center. Ra is Right ascension and Dec is Declination in 
the equatorial coordinate system. Since the data set is within 
2 < 1, the radial component can be ignored, and Dc{z) is 
written as 


D, 

E{z) 


c r dz' 

'id Jo ^ 

'\/ Dm(l “1“ 2 )^ -f Da 


( 2 ) 

( 3 ) 


Here, we assume the universe is flat. flA is the dark 
energy component and Qrn is the matter component. Da + 
Qrn = 1. The catalogue includes two sub-catalogues (co¬ 
moving and redshift) according to the coordinate in which 
it is searched for the voids. Voids in comoving catalogue 
are found out through assuming dark matter component 
Da = 0.73 and Voids in redshift catalogue is found out 
through assuming dark matter component Da = 1.. The 
void number and redshift range in each section are listed in 
table [T] We use only the ’’central” void s which are selected 
to avoid the boundary or mask effect IlSutter et al.M2012bl . 

mm- 


3 TESTING ISOTROPY AT LARGE SGALES 

The void alignment depends on three sources: initial con - 
ditions, void packing and tidal forces dPlaten et ^ I 2 OO 8 I I. 
These three sources should be local effects and should not 
contribute much at large scales. Therefore the distribution 
of the orientation of the cosmological voids should be a ran¬ 
dom distribution. Based on this, it was suggested that the 
void orientation with respect to the line of sight could be 
a measure of the reds hift distortion llFoster fc Nelsor]|2009l : 
iRvden fc Melottlll996h . 

If void orientation is random, the angle, 9, between its 
orientation direction and line-of-sight satisfies the distribu¬ 
tion 

p{e)de = smede (4) 


or 


P{cos9)dcos9 = dcosO 


( 5 ) 


To examine this assumption, we calculate the orienta¬ 
tion of voids through ellipsoid approximation. We compute 
the location of the void’s center through the volume weight 
as 




( 6 ) 


Here, T* is the position of the void’s i-th galaxy’s and 
Vi is its Voronoi cell’s volume. We then compute the shape 
tensor, Sij, 


Sij = '^{Xi - Xci){xj - x^j)mk 

k 


( 7 ) 


k is the void galaxy’s index, i and j are the three po¬ 
sition’s indexes, mk is the mass weight. Here we choose 
ruk = 1 for all the cases. The three eigenvectors are 
the orientation of the ellipsoid’s axes. The three eigen¬ 
values are their axes’ length square. We calculate the 
orientations under Da = 0.73 to compare with recent 
WMAP and Planck’s measurem ent dKomatsu et al.l 1201 ll : 
IPlanck Collaboration et al.ll20ljl . We use voids in the co¬ 
moving catalogue and show the result in figure [T] At first 
glimpse, one finds that the long axes have values larger than 
1 for small cos 6. That means the long axes are anti-aligned 
to the line-of-sight. It is not very clear whether median and 
short axes have any preferred direction without a detailed 
analysis. 

There are at least two reasons to this a nti-alignment 
phen omenon. One is Alcock-Paczynski effect jAlcock et aP 
I 1979 II and the other is peculiar velocity. We will discuss them 
in the next section. 


4 TESTING ALGOGK-PAGZYNSKI EFFECT 
ON VOIDS’ LONG AXIS AND PECULIAR 
VELOCITY AS A CAUSE OF STRETCHING 

It is known that Alcock-Paczynski effect could cause redshift 
distortion. At the same time it also distorts the orientation of 
a void. For our purpose, to study the Alcock-Paczynski effect 
on the orientation, we can’t approximate the void shape with 
a sphere, because of the voids’ irregular shapes and their 
tracers’ peculiar velocities. To simplify the argument, we 
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Figure 1. The three axes’ orientation distribution. We set = 
0.73 for the figure. We use voids in the conioving catalogue in 
the 9 sections. The orientation of the long axes are slightly anti¬ 
aligned to the line-of-sight. The bin size of cos6 is 1/30. 


consider the void to be a one dimensional rod. The rods’ 
directions are the orientation of long axes. The real situation 
is more complicated, but our approximation can at least tell 
us what is going on when considering the Alcock-Paczyski 
effect on the void’s orientation. 

Let’s consider a one dimensional void. Its angular extent 
is S9v and redshift extent is Sz^. In the comoving coordinate 
system, 


drc = 59vDc ( 8 ) 

5dc = SzvdzDc (9) 


Sdc void’s comoving size in line-of-sight direction and 
Svc is the void comoving size in angular direction. Therefore, 
the angle between the void’s long axis and line-of-sight is 


tan 0c 


Dc S9y 
dzDc dzv 


( 10 ) 


0 c follows a random distribution if Dc is the real co¬ 
moving distance. If one uses the other factor, D„, {Dn = 
in general discussion), one will have mistaken the angle be¬ 
tween the void’s orientation and line-of-sight as 


tan0„ 


Dn 59 n 

tdz Dn 5Zv 


( 11 ) 


In the measurement, Szv and 59^ are measured. How¬ 
ever, the comoving distance, Dc, is not known directly. One 
may use instead of Dc to study the voids’ orientation. 
The relation between 9n and 0c is 


tan0„ = /( 2 )tan 0 o 

(12) 

. dzlnDc 

fV) “ o , n 

Oz In Dn 

(13) 

The distribution 



(14) 


In the actual isotropic coordinate, the void’s long axis’s 
are random distributions and therefore P{9c) = sin 0c. From 
equation m one finds 


Figure 2. cos{9n) distribution. It aligns to the line of sight as 
/>!• 


P{9n) = 


sin 9n 
~1^ 


^ cos^ 0„ -b 


sin^ 9n \ 

~1^ ) 


3 

2 


(15) 


This is the angular distribution function for a coordi¬ 
nate different from the comoving coordinate system. From 
figure m one finds it is anti-aligned to the line of sight as 
/ > 1 and aligned to the line of sight as / < 1. 

Although the voids’ shapes are complicated, we could 
still find that the voids’ long axis share the same effect. To 
examine the effect, we use drlOcmassS dataset because of 
its higher redshift, and vary Qa from 0.0 to 1.0. Figure [3] 
shows that the voids’ orientation changes from anti-aligning 
to aligning to the line-of-sight as we go from Qa = 0 to 
SIa = 1 . Based on this, one expects that one could find the 
dark energy and dark matter component (Ha and flm) by 
searching for P{9n) which best satisfies a random distribu¬ 
tion. For a pure random distribution P{9), < cos0 >= 0.5. 
Therefore, one could calculate the average < cos 9 > to 
check whether the distribution follows a random distri¬ 
bution. Figure |4] shows how < cos0 > changes with re¬ 
spect to Ha. The figure shows that Ha = 0.78 ± 0.02 in 
la interval. Howeve r, we must point out that an empiri¬ 
cal st retching factor llSutter et aHl2014a : iLavaux fc Wandeltl 
l2012tl . which may be caused by the peculiar velocity, must 
be in cluded. The same effec t has been noticed in the SDSS 
DRh dFoster &: NelsonI [ioO^ L They found that the average 
orientation angle increases to higher 0 values as z from 0.1 
to 0.16. If the empirical stretching factor doe s exist, then 
1/f = 1.16 ± 0.04 in a comov ing coordinate dSutter et al.l 
l2014d : iLavaux fc Wandeltll2012l ~) and it will not be isotropic 
even in the correct comoving coordinate. We can test the em¬ 
pirical stretch factor by calculating the average < cos0 >. 


cn/1 

<cos0>= / P(0)cos0d0 (16) 

Jo 

For 1// = 1.16, < COS0 >= 0.46. We plot all < cos0 > 
in Ha = 0.73 case. Figures [5] shows that most of < cos0 > 
are less than 0.5, which is the value of a pure random distri¬ 
bution. This confirms the existence of the extra stretching 
factor. This extra stretch is from the void galaxies’ pecu - 
liar velocities dSutter et al.|[2014d : ILavaux fc Wandeltll2012l L 
However, the redshift type of a void and comoving type of a 
void do not give exactly the same result. This confirms that 
the void galaxies’ quantities highly depend on the coordi- 
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Figure 3. The probability function changes with different Qa- 
Here we use drlOcmassS dataset. 




Figure 4. < cos0 > in different The black line is the ex¬ 
pected < cos 0 > in the isotropic case. The vertical dash lines are 
the 1(7 interval of possible Qa- The data point are redshift types 
in drlOcmassS section. 


nate which is used in the void finder (|Nadathur &; Hotchkiss! 
|20 i4) . 


5 VOIDS’ UNIVERSAL DENSITY AND 
REDSHIFT DISTORTION 

The previous section has confirmed that the voids’ orienta¬ 
tion depends not only on redshift effect, bnt also on an em¬ 
pirical factor. According to simulations, the em pirical factor 
is cau s ed by the galaxies’s pec nliar velocities lISutter et al.l 
l2014d : iLavaux fc Wand^l2012l ~). In general, it is very dif¬ 
ficult to study peculiar velocity without any other assump¬ 
tion. Howeve r, it has been found that voids’ density profile 
is universal ( Hamans_et_aJj 20lJ; |Layau2Lj£_Wandehj 20l3; 
Colberg et alj 20051 : Padilla et al.l 20051 : Ricciardelli et all 


20131 . 20 14) and its average v elocity fit the l inear theory 


very well I Hamaus et al.l l2014l : |Paz et ahl l2013ll . Therefore, 
we could study the peculiar velocity from this universal den¬ 
sity profile. 

Here we adopt the density profile form from paper 


Figure 5. < cos0 > in different catalogue section. The black one 
is from the redshift type of the catalogue, and the red one is from 
the comoving type of the catalogue. f^A = 0.73 in the calculation. 
The black line represents the random distribution result and the 
dash line represents a stretch factor 1// = 1.16. 


(|Hamaus et al.|[2014h . 


Pvir) 


— 1 + 


1 - {r/vsY 
1 -I- {r/r^y 


--nZ^HvAir) 


AW = 


(17) 

(18) 
(19) 


Here, 7 = 0.55. r = \/x?+Xo + x^. is velocity 


from the linear theory dPeebled 197(il i. It has been found 
that a a nd (j a,t z = 0 are related to rs/r„ in the follow¬ 
ing form (|Hamaus et al.ll20lA . 


P = 


= -2{rslr^ - 2 ) 

17.5rs/r„ — 6.5 for rs/r„ < 0.91 
—9.8rs/r„ -|- 18.4 for rs/r„ > 0.91 


( 20 ) 

( 21 ) 


However, we study voids in z y 0 case and therefore 
we should treat them as two free parameters. We hope that 
a and P dependence on the redshift could be found in the 
future studies. So far, *3 can not be obtained directly from 
redshift, and one must consider the distortion effect from 
the galaxies’ peculiar velocities. Therefore, 


X3 = Z — v/H 
c 


Z = 


H 


-5z 


( 22 ) 

(23) 


where v is the galaxy’s velocity and 5z is the redshift differ¬ 
ent between galaxy and void center. Since the galaxy’s ve¬ 
locity is involved in the measurement, one needs a velocity 
marginal function to find x^. The simplest marginal function 
will be a Gaussian distribution with the center at Vv How¬ 
ever, since we rescale the void according to its effect radius, 
{Rv = , V is the void’s volume) while stacking voids, 

the velocity will be rescaled too. Therefore whether Gaus- 
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sian distribution is a velocity distribution function is un¬ 
clear. We hope future N-body simulation could also provide 
the form of marginal function. Right now we use Gaussian 
distribution. 


fgiv) 


1 

y/^Vo 


exp(- 



(24) 


Vo is the root mean square of the velocity. Exponential 
distribution is also a popular distribution in two-body cor¬ 
relation. We had tested it and found the result is similar 
to Gaussian distribution one, therefore we are going to give 
only Gaussian distribution’s result. Of course vq should be 
a r dependent function. Finding its actual form will rely on 
the future studies, and here we will treat it as a constant. 
The average density in {xi, X 2 , Z) is 


3.2x10"'- 

■ ■ 

2.8x10"'- ■■ 

■ 

■ ■ ■ 

,2.4x10"'- . 

lo." 

■ ■ 

2.0x10"'-* ■, 

■ ■ 

1.6x10"'- 

0.48 0.52 0.56 0.60 


Po{xi,X2,Z-nA)dZ = f 

J P r 

S{Z — X3 ——l^)dx3dw (25) 

where 5{ri) is a Dirac Delta function. Da appears in the 
formula, because Hubble constant and position are functions 
of Da. Different Das give different mass density. We first 
rescale the voids given in comoving coordinates according to 
their effect radius (7?„) and then stack them together. We 
expect a single galaxy’s peculiar velocity could be several 
hundred to thousand km/s. This will give several to dozens 
Mpc of distortion in line-of-sight (Z) direction. Therefore we 
focus only on voids bigger than lOMpc. Low redshift voids 
either do not satisfy or the void number is very low (less 
than one hundred). Considering the void number, we choose 
drl0mass2 to study the effect. We use voids in comoving 
coordinates to minimize the void finder’s selection effect. 

In order to find the density profile, one needs the av¬ 
erage galaxy density. This however is not easy, because the 
average galaxy density depends on survey’s capability and 
is redshift and location dependent. Figure [S] shows how the 
galaxy density changes with respect to redshift in the co¬ 
moving space. It is clear the galaxy density highly depends 
on the redshift. At this moment we tr eat this density as aver¬ 
age density (Pm(z)) at eac h redshift llNadathur fc Hotchkiss! 
l2014l : l^adathur et al.ll20l3l . We then stack the voids accord¬ 
ing to void radius, Rv The density profile is calculated in 
the following way 


pa{Z^ dv] Da) 


Pi(-^, d„; Da) 

W pmlZi) 

i 


(26) 


Pi{Z, Da) is the galaxy density of ith void at {Z, d„) 
re spected to Da. Here Zi is the redshift of the bin. dv = 
y/Xi + X 2 and N is the number of voids in the stacking. 
The error of each bin is calculated from 


a{Z, Da) = 


1 


N{N-1) 




PilZ^dv\^A^ - ,2 

— Pa) 


pm Vi} 


(27) 


Figure [7] shows the stacking density profile according 
to different void size. We stack 175 voids for > 40Mpc 
case, 202 voids for 40Mpc > Rv > SOMpc case and 236 
voids for SOMpc > > 20Mpc case. The bin size is chosen 

to be 3Mpc X 3Mpc. The stacking void number does not 
represent the stacking galaxy number in study because the 


z 

Figure 6. The galaxy density according with respect to redshift: 
The density is calculated according to Da = 0.73. The density 
depends highly on redshift and can not be treated as a constant. 
The density unit is Mpc“® 


2.0-1 



0.0-1—^^.-,- ^^^^^^-, 

0 20 40 60 80 100 

R 


Figure 7. The voids are separated into three different radius 
range (20Mpc < Rv < SOMpc, SOMpc < Rv < 40Mpc and 
40Mpc < Rv) stacking according to their radius. We rescale all 
the voids to SOMpc and stack them together. The bin size is 
SMpc X SMpc in {Z,dv) coordinate. 


different void has different volume and there are different 
masked regions to avoid. The smaller voids’ profile has a 
higher density right outside the voids’ radius. Apparently 
the density profile depends on the voids’ sizes. Therefore it 
is unlikely to stack each void without considering its radius. 
At the following we study void density profile according to 
these three stacking void profiles. Also we choose 5^ = — 1, 
because there is almost no galaxy near the void center. 

Figure |8] shows average density(pa(^, du)) and error 
{a{Z,dv)). Near Rv = 40Mpc, there is a point with much 
lower pa{Z,dv) and its a{Z,dv) is much less than its neigh¬ 
borhood. This shows that if the average density {pa{Z, dv)) 
is much smaller than average galaxy density, then a{Z,dv) 
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3.0 


20<R <30 



100 



Figure 8. The density distribution with error bar: The stacking 
voids’ effect radius are from 20Mpc to 30Mpc. = 0.5. Since the 
limitation of dataset, sometimes the error bar is underestimated. 
In this case, a point near Rv = 40Mpc is much smaller than its 
neighborhood. This causes that its error bar is underestimated. 



Figure 9. The drl0mass2 density distribution: we stack 400 
voids. The voids’ effect radius are from 30Mpc to 40Mpc. The 
bin size is 3 x 3Mpc^ in {Z,dy) coordinate. Here Qa = 0.73. We 
choose Re = 30Mpc in the figure. 


is underestimated. If these underestimated points are in¬ 
cluded in the calculation the increase in will be signif¬ 
icant. Since these abnormal points appear at small dy and 
are caused by the limitation of the dataset, we therefore re¬ 
move dv < 3Mpc data points to avoid the underestimate of 
(7{Z, d„). 

Figure [5] shows the drl0mass2’s average density profile 
in Z — dv plane. (Here we show = 0.73) One finds exten¬ 
sion in Z-direction at small radius(~ lOMpc) and stretch at 
larger radius (~ 40Mpc). Highest density peak appears near 
the Z-axis. This is an evidence of galaxy’s peculiar velocity 
causing the density distribution distortion. 

We then use the universal density profile to fit the stack¬ 
ing data. We assume each bin is independent from the oth¬ 
ers. Their cross correlations are independent and diagonal. 


Figure 10. The density profile by fitting empirical density func¬ 
tion. The other parameters are the same as figure [9] 


Therefore we include only ij^{Z, IIa) to calculate x^- 'C 


is calculated according to 




{pa{Z, dy) - Po{Z, dv)Y 
(j'^{Z,dv-,Q.A) 


(28) 


We use data between lOMpc to 60Mpc to calculate 
Figure [TU] is the fit of the density profile from the empirical 
density profile. The result shows similar density stretching, 
extension and distortion as in figure |9l 

One of our goals is to find the cosmological dark energy 
component, Ha. Figure ITT] shows the minimum x^ changes 
with respect to Ha (The other parameters are not fixed). We 
fit 291 points with 5 parameters. Therefore a reasonable x? 
value is ~ 286 ± V286. x^ at small Q,a is farther away from 
this region and can be ruled out. In 30Mpc < Rv < 40Mpc 
and 40Mpc < Rv case, the x^ minimums is still a little too 
large and their best fit parameters are not reasonable (uo less 
than lOOkm/s which is unreasonable because the bin size, 
3Mpc, at least induces 300km/s uncertainty on velocities). 
This implies a further study on the form of universal density 
profile is needed. At the same time, we assume all the points 
are independent and the likelihood is 


L{Qa, a, /3, Wo, rs) = b[ 


(27rcr2)l/2 


exp(-xV2) 


(29) 


The probability for a particular Qa is calculated accord¬ 
ing to the following marginalization 

P{Qa) = j L{QA,a,h,vo,rs)M{flA,a,fi,vo,rs)dadpdvodra{30) 

Here M is a prior probability function. It depends on 
several parameters, but in general M is assumed to be a 
constant. Figure [T^] shows the likelihood ratio under this as¬ 
sumption. Again the small Has are ruled out immediately. 

But we must emphasize that the minimum x^s appearing in 
30Mpc < Rv < 40Mpc and 40Mpc < Rv cases are not rea¬ 
sonable and therefore M cannot be taken as only a constant. 

A further study of the prior probability is needed. 


© RAS, MNRAS 000,[iH8] 






































Void alignment and density profile applied to measuring cosmological parameters 7 



Figure 11. The vs. Here 13, vq and rs are not fixed. 


0 - 

- 20 - 

-40- 

-o -60- 
& 

Si 

1 -80- 
- 100 - 
-120 
-140 


vr 


/ 


20<R/30 

30<R^<40; 

40<R, 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


caused by the peculiar velocity, it is more realistic to apply a 
model which can quantify the peculiar velocity effect. It has 
been pointed out by Hamaus etal. that voids’ density profiles 
are universal both inside and outside voids and their aver¬ 
age v elocities satisfy linear models very well dHamaus et ^ 
l2014li . Therefore we apply Hamaus’s density profile to the 
voids density profile in SDSS lODR. 

Although it is ultimately better to stack all the avail¬ 
able data at once, figure [7] shows that the density profile 
also depends on the void’s effective radius. This reduces the 
amount of voids that can be stacked together at once. In 
our study we separates the voids into three different effec¬ 
tive radius intervals and fit the void density profiles with the 
empirical density distribution form. The result is consistent 
with the tsst and it implies that the empirical density 
profile form we obtained is reasonable. However, this is not 
enough to apply the result to calculate the confidence level 
directly, since the prior functions must be known. Unfortu¬ 
nately, in this case assuming a constant prior function is not 
suitable. It is unclear if this is caused by the limitation if 
the voids’ number or by the voids’ density profiles. This in¬ 
dicate that further analysis is warranted in order to better 
understand the voids’ statistics and associated phenomena. 
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Figure 12. The probability ration vs. We marginalize a, /3, 
Vq and Vs parameters, and assume the prior probability function 
is a constant. Pq is the maximum probability. 

6 CONCLUSION 

In this paper, we studied the properties (orientation and 
density profiles) of the cosmological voids using the SDSSIO 
data. We first confirmed that the voids’ orientation is not 
purely random as one might naiv ely expect. This effec t 
has been first noticed by Foster etc JFoster &: Nelsonll2009li . 
Later, Sutter etc. have found a constant empirical stretch¬ 
ing factor using numerical simulation s, and applied the 
stretching factor to SDSS DRIO data dSutter et aP l2014cl : 
iLavaux fc Wandeltll201 j l. They found that the voids’ ori¬ 
entation can be used as a dark energy probe. Along these 
lines, a constant stretching factor will make the study of 
Alcick-Paczynski effect relatively straightforward. However, 
as illustrated in figure [4] we showed that the stretching fac¬ 
tor highly depends on the coordinate w hich is used in the 
void finder dNadathur fc Hotchkissl20i3 l. This is not surpris¬ 
ing, because different void finders give different void galax¬ 
ies. Because of this discrepancy, it is very difficult to study 
the Alcick-Paczynski effect directly through the void orien¬ 
tation. To reduce this discrepancy, we included the galaxies 
outside voids into our analysis. Since the stretching effect is 
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